Robustness of d-Density Wave Order to Nonmagnetic Impurities 
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Effect of finite density of nonmagnetic impurities on a coexisting phase of d-density wave (DDW) 
order and d-wave superconducting (DSC) order is studied using Bogoliubov-de Gennes (BdG) 
method. The spatial variation of the inhomogeneous DDW order due to impurities has a strong 
correlation with that of density, which is very different from that of DSC order. The length scale 
associated with DDW is found to be of the order of a lattice spacing. The nontrivial inhomogeneities 
are shown to make DDW order much more robust to the impurities, while DSC order becomes very 
sensitive to them. The effect of disorder on the density of states is also discussed. 
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a. Introduction One of the recent proposals in the 
context of high temperature cuprates is that a true 
broken symmetry state dubbed as d-density wave state 
(DDW) is responsible for the pseudogap phenomena. |l| 
This phase was first suggested in relation to the exci- 
tonic insulators 2] , and it was found as one of the ground 
states of the t-J type model. 3] The DDW is a parti- 
cle hole condensate with angular momentum 2. The or- 
dered state can be characterized by the circulating cur- 
rent arranged in an alternating pattern on a square lat- 
tice, which can be detected as a Bragg scattering signal in 
neutron scattering measurements. [3, |5j But the neutron 
scattering experiments 0, 0, 13 in cuprates remain con- 
troversial. Thus, the definite conclusion on the relevance 
of the DDW order to the cuprates requires more precise 
experiments on various doping concentration of cuprates, 
and further theoretical studies on the properties of this 
new order. Especially, the effect of the nonmagnetic im- 
purities on DDW order is an important subject to inves- 
tigate, since any well-prepared cuprate sample contain 
an intrinsic disorder, minimally from non-stoichiometry. 

The simplest possible description of the impurity effect 
is the self-consistent T-matrix approximation (SCTMA) 
0. This mean field picture excludes not only the free- 
dom of the ordered patterns, but also the interference 
of the impurities. Within this approximation, the ther- 
modynamics were found to be identical to those of a d- 
wave BCS superconductor (DSC) in the unitary limit. fl0| 
From the density of states, one can see that electrons are 
localized close to the Fermi energy, and the change in the 
transition temperature is given by the Abrikosov-Gorkov 
formula known in BCS superconductors. [lOj Within the 
standard non-crossing approximation, the similarity be- 
tween the DDW and DSC is based on the d-wave sym- 
metry of the gap. 

In this paper, we study the effect of impurities on 
DDW order and for the case where DDW coexists with 
DSC using Bogoliubov-de Gennes (BdG) technique. This 
method is the mean field approximation, but it allows 
spatial inhomogeneity in order parameter. In the case of 
the disordered DSC with a short coherence length, it was 
shown that the superfiuid stiffness is significantly larger 



than that predicted by the SCTMA, due to the nontrivial 
spatial structures of the order parameter|ll). 

We found that the DDW order is more robust than the 
DSC order to the impurities, which cannot be understood 
within the conventional T-matrix approach. The physical 
ground for our findings will be discussed later. 

b. Model We model two dimensional disordered 
DSC and DDW order by the following Hamiltonian. 

(i) 

<ij> <ij>,a,f3 

The first term is the kinetic energy which describes elec- 
trons, with spin a at site i created by c\ a , hopping be- 
tween nearest-neighbors < ij > on a square lattice. The 
disorder potential V(i) in the second term is an indepen- 
dent random variable at each site which is either +Vq, 
with a probability nj mp (impurity concentration), or zero, 
and \x is the chemical potential. The last, interaction 
term 12] is chosen to lead to a coexisting DSC and DDW 
order ground state in the disorder-free system, where Si 
and rij are the spin and density operators, respectively. 

The mean field decomposition of the above Hamilto- 
nian leads to following BdG equations [lH Il5| . 

(1. < 2 » 

where £u n (j) = -J2 s {t + *(.?; 6)e~ i( *- r i }u n (j + 5) + 
{V{j) r fij)u n {j) and Au n (j) = £ tf A(j + 8;5)u n (j + 6), 
and similarly for v n (j). The local DSC pairing (A) and 
DDW (x = Im^f) amplitudes on a bond (j; 6) are defined 
by 

A / • J +W . 

A u;<>) = 1 — ( c j+sicn + c 3 i c j+si), 

*(r,s) = J - J ^{c] +Sa c Ja -c J J+Sa )e-^Ai) 

where 6 = ±x, ±y. The inhomogeneous Hartee and Fock 
shifts are given by p,j = \x + (j + W) J2s( n j+s) an( ^ 
i?e[^(j;5)] respectively. 



2 




\ 

\ 

\ 

\ 
\ 



\ (e)—~ a— A/t (with finite x. \ 

\ (n>=0.95) 

Q L I I I I 1 I I I I i I I i 1 i i 1 1 — ^ — i 1 > ' i 

0.01 0.02 0.03 0.04 0.05 

FIG. 1: The evolution of A and x (normalized by their pure 
values) with n llnp is plotted under different conditions (see 
text). 



We numerically solve for the BdG eigenvalues 
E n > and eigenvectors (w n ,t>„) on a lattice of N 
sites with periodic boundary conditions. We then 
calculate the d-wave pairing amplitude A(j;S) ~ 
(J + W)J2 n K(3 + S)v* n (j)+u n (j)v*U + S)]/4: 
and the DDW order and Fock shift as the 
imaginary and real parts of ^(j; S) = ( J + 
2W) E„ KUM +8)- MjXU + S)] /4 at T = 0, 
and the density (rij) — 2J2 n \ v n(j)\ 2 - These are fed back 
into the Bd G eq uation, and the process iterated until self 
consistency [l6| is achieved for each of the (local) vari- 
ables defined on the sites and bonds of the lattice. The 
chemical potential /i is chosen to obtain a given average 
density (n) — J2i( n i)/N- We define the site dependent 
order parameters in terms of the bond variables as, 
A(j) = [A(j; +4) - A(j; +y) + A(j; -x) - A(j; -y)] /4 
and similarly for x(j)- 

We have studied the model at T = for a range of pa- 
rameters and lattice sizes up to 40 x 40. Here we focus on 
J = 1.16, and W = 0.6, in units of t = 1, with (n) = 0.95 
on systems of typical size 30 x 30. For these parameters, 
and rii m p = 0, the maximum DSC gap is A max = 0.16 
and the maximum DDW gap is Xmax = 0.31. In the 
pure system our calculations reproduce a phase diagram 
of A max and Xmax as functions of filling similar to Ref. 
[l3j . For the impurity potential we choose Vq = 100, 
close to the unitary limit. The results are averaged over 
10 different realizations of the random potential. 

c. Effect of Impurity on DDW and DSC Orders We 
summarize our main results in Fig. (1), where we plot 
the disorder dependence of different orders (normalized 
to n; mp = values). Let us first look at the line (a) 
that represents the behavior of x as a function of rii mp at 



half filling ((n) = 1). At this filling DDW is the stable 
order and DSC order in fact vanishes for the pure sys- 
tem. Comparing the behavior of \ with the results from 
SCTMA calculations (represented by (c) curve) we see 
that the DDW order is more robust to impurities than 
predicted by SCTMA. 

On the other hand, away from half filling when we 
force x = in BdG equations, DSC becomes the sur- 
viving order and the n; mp dependence of A is given by 
the (b) line, which in fact is very similar to \ in the (a) 
curve. Such robustness of the DSC order to impurity had 
been studied before 0| , and it is attributed primarily to 
the fact that - each impurity affects superconductivity 
rather inhomogeneously by destroying SC order within 
a small region (of size determined by coherence length 
£) around it. Hence the long range order is not glob- 
ally affected. We find from our current numerical results 
that similar picture holds for DDW order as well, and x 
is also affected locally by impurity, keeping long range 
DDW order robust. 

However, similar study for the coexisting phase of DSC 
+ DDW order at (n) = 0.95 reveals surprisingly that, 
superconducting order is severely affected by disorder 
(curve (e)) in the coexisting phase, much more so than in 
the absence of x- On the contrary, the DDW order (curve 
(d)) coexisting with DSC order becomes even more ro- 
bust to impurities. In fact for low n lmp , x even increases 
with impurity. The rest of the paper is organized towards 
the detailed understanding of these unexpected results. 

From Fig. (Id) we saw that the DDW order increases 
for small n; mp . To get a further insight, we study the 
spatial structures of the order parameter (particularly at 
large n; mp ) on the lattice for each impurity configuration. 
In Fig. (2a) we present a Grey-scale plot of x on a typical 
30 x 30 lattice at rti mp = 0.06 for a given realization 
of scatterers. The dark (light) regions represent larger 
(smaller) values of x- Comparing this structure with Fig. 
(2b), that gives the spatial distribution of \(n) — 1| for 
the same ni mp , we see that x is large in space where 
local density is close to 1 (half filling). The strong spatial 
correlation between these two panels is striking, although 
it is not exact; the scale of modulation of x is somewhat 
larger than that of density. However, the strong tie of 
local (n) and x suggest that the length scale of fluctuation 
of x would be governed by that of (n), which is rather 
small (of the order of hp 1 ). This can be understood as 
follows. 

The length scale associated with the DDW order, 
£ddw ~ 1/x- When impurity is introduced, the bond 
current attached to the impurity site is forced to be zero. 
So does the density. However, the bond current "near" 
the impurity site is re-constructed to satisfy the current 
conservation, and one should note that the healing length 
is of order of a lattice spacing. How the magnitude of the 
re-constructed bond-current is determined? This magni- 
tude is strongly related to the local density. The electron 
density depletes close to impurities and increases at lo- 
cations far from it, to keep the average at the desired 
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FIG. 2: Grey scale plot of DDW and DSC order on lattice 
for riimp = 0.06 and for a particular configuration of unitary 
impurities. The dark (light) region indicates large (small) 
values of the variables on given locations. Panel (a) and (b) 
are the \ and | (n) — 1 1 respectively for a system with coexisting 
DDW + DSC order at (n) = 0.95. A is shown in (c) with a 
finite x, and (d) with \ forced to zero everywhere. 



value. Since at low disorder, a large number of sites at- 
tain (iii) ~ 1, x increases at those sites; the DDW order is 
most stable near half filling, where perfect nesting occurs 
for our model. As a result average \ increases. At very 
large n; mp , local density would be either much larger or 
smaller than 1, and \ would decrease everywhere. This 
argument can be substantiated by looking into our results 
for each configuration of impurities. 

For (n) = 1 , introduction of impurity makes local den- 
sity only to deviate from half filling. As a result \ de- 
creases monotonically as found in Fig. (la). The above 
argument for the behavior of DDW order with impurity 
is independent of the coexisting DSC order and we also 
found similar trend in \ a s in Fig. (Id) for (n) < 1 even in 
the absence of DSC order, which is consistent with our 
picture. This shows that DDW order responds to the 
density fluctuations due to impurities. On the contrary, 
the DSC order in the presence of impurities is not related 
to the local density fluctuations as DDW is, even though 
the length scale, £dsc ~ which is of the order of a 
few lattice spacing for high temperature superconductors 
under considerations. The behavior of A in the presence 
of impurities is shown to be related to the electron-hole 
mixing in the real space |l5l | ; A is large when the density 
is close to the chemical potential. 

Fig. (2c) and (2d) presents the spatial structure of A 
on lattice with the same rii mp configuration in the pres- 
ence and absence of DDW order. We clearly see that the 
DSC is strongly suppressed by the impurities when coex- 




FIG. 3: Panel (a): Disorder averaged density of states N(u) 
on a system with only DDW order (at (n) = 1) on 10 x 10 
unit cells each of which are of size N = 30 x 30. It shows 
the robustness of DDW order (see text). Panel (b): Similar 
to Panel (a), but with coexisting DDW + DSC phases (at 
(n) = 0.95). With DSC gap is washed out, 

however, strong DDW order persists. 

isting with DDW order (as also observed in Fig. (lb) and 
(le)). The existence of DDW strongly affect the strength 
of the DSC, because away from the impurities there are 
regions where the density is near half-filling, hence the 
DDW becomes strong. Strong DDW allows significant 
weight of (7r,7r) scattering that mixes the + and — lobes 
of the DSC order and thereby DSC becomes weak. The 
regions of small density does not contribute to DSC or- 
der as well, due to the absence of enough electrons for 
pairing! Thus in the inhomogeneous coexistence phase 
DSC order is suppressed everywhere. 

d. Averaged Density of States Let us now study 
the (impurity) averaged density of states (DOS) 
= Jr E„,i [\u n (i)\ 2 S(u; - E n ) + \v n (i)\H(u + E n )} 
(where we broaden the delta functions with a width com- 
parable to average level spacing). To obtain a better 
statistics for N(uj), we used "Repeated Zone Scheme", 
(Tt| that describes a large effective system made out of 
10 x 10 unit cells, each of which is of dimension 30 x 30. 
In Fig. (3) we plot N(u>) as a function of oj for differ- 
ent rtj mp for the case of only DDW order (Panel a) and 
coexisting DDW + DSC order (Panel b). For the pure 
system with only DDW order, N(u) is the standard d- 
wave DOS. With increasing nj mp we see that the gap-edge 
singularities get rounded off and a small accumulation of 
states is produced at the particle side of spectrum close 
to lo = 0. The accumulation of electrons around a sin- 
gle impurity effectively provide impurity screening, which 
will produce enhanced states at the particle side of the 
spectrum. |l8j| Such resonances from each impurity con- 
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tribute to the average N(u>) and produce a broad band 
which is reflected as a bump in Fig. (3a). However, the 
strength of the DDW order is not affected much (given 
by the relative location of the two coherence peaks). At 
this point, we should emphasize that the DOS structure 
for impure DSC state is very different, where coherence 
peaks get strongly suppressed and a thin gap persists at 
v = [H Il9|. so that N(0) = for all n imp . From 
our results with DDW order, we find that iV(0) oc ni mp , 
which is in disagreement with the prediction of T-Matrix 
result (N(0) <x ^n^MM- 

In Fig. (3b), for coexisting DSC + DDW, a double-gap 
DOS is expected at rt imp = [2(J; superconducting gap at 
lo = and d-density wave gap at u = fi. With increasing 
"-imp, DSC gap gets washed out and by n lmp = 0.03, N(u>) 
looks very similar for Fig. (3a) and (3b) (The overall shift 
for the later case is due to the particle-hole asymmetry). 
This demonstrates in a different way our main result, 
that, the DSC order is very sensitive to impurity whereas 
DDW order is robust in the coexisting phase. 

e. Summary and Discussion We studied the effect 
of nonmagnetic impurity on DDW ordered state using 
BdG technique. While the standard SCTMA indicates 
that the effect of impurity on DDW is similar to that on 
DSC, we found that the spatial variation of the DDW 
order has a strong correlation with that of density [Fig. 
(2a) and (2b)] , and it is very different from that of DSC 
order [Fig. (2a) and (2c)]. We discussed that this occurs 
because the length scale associated with the DDW order 
is of order of a lattice spacing (~ I/Hf), which suggests 
that the spatial variation of DDW order is related to the 
density fluctuation, while the DSC order is related to 



particle-hole mixing. Therefore, the effect of impurity on 
the DDW order is very different from that of DSC order, 
which can not be obtained from the standard SCTMA 
method. 

When DSC and DDW coexist, it turns out that DDW 
order do not care about the existence of DSC and it still 
follows the density profile in the presence of impurity. 
However, DSC order would vanish almost everywhere 
[See Fig. (2c)]. This is because in the region of larger 
density it is killed by DDW, and in the region of smaller 
density it is destroyed by disorder. Thus in the inho- 
mogeneous media both DDW and impurity are acting to 
suppress the DSC order. 

Our current picture brings out the unexpected results 
and their understanding at the mean field level; if the 
DDW phase exists in cuprates, the Bragg signal would 
be detected in neutron scattering measurements even in 
the presence of strong nonmagnetic impurity, while the 
width of the Bragg peaks depends on strength of impu- 
rity. However, the definite answer for its relevance to the 
cuprates requires the understanding of the role of strong 
correlation, and interplay between different competing 
orders, which warrants further studies. 
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